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(I)  INTRODUCTION 


This  final  technical  report  is  concerned  with  the  tasks  specified 
in  the  governing  contract.  They  are  as  follows: 

(1)  Computer  solutions  developed  during  the  previous  contract 
year  will  be  used  in  making  a short  (one  to  two  minutes) 
movie  film  Illustrating  unsteady  flow  with  a shock  wave  in 
a two  dimensional  symmetric  channel;  the  unsteady  behavior 
will  arise  as  a result  of  oscillatory  pressure  pulses  intro- 
duced in  a plenum  downstream  of  the  shock  wave. 

(2)  The  study  of  the  characteristic  time  for  the  formation  and  col- 
lapse of  a shock  induced  separation  bubble  will  be  completed. 

(3)  An  attempt  will  be  made  to  use  methods  developed  previously 
to  develop  a criterion  for  incipient  separation  Induced  by 
an  oblique  shock  wave  impinging  on  a flate  plate  turbulent 
boundary  layer  in  transonic  and  in  supersonic  flow;  the  ex- 
tent to  which  this  analysis  will  be  carried  out  will  depend 
upon  the  time  remaining  after  completion  of  the  first  two 
(2)  tasks. 

The  solutions  referred  to  in  task  (1)  are  shown  in  reference  1. 

They  describe  the  motion  of  a shock  wave  and  the  time  and  space  varia- 
tions of  the  flow  properties  downstream  of  the  shock  wave  in  a nozzle 
flow  when  the  pressure  in  the  plenum  (back  pressure)  is  made  to  oscil- 
late. One  of  the  more  Interesting  cases  for  which  they  are  valid  is 
that  in  which  the  shock  wave  moves  up  to  and  passes  through  the  nozzle 
throat,  disappearing  upstream  in  the  subsonic  flow.  The  solution  for 


the  shock  wave  velocity  given  in  reference  1,  being  a first  approxima- 
tion, gives  an  infinite  velocity  for  the  shock  wave  as  it  passes 
through  the  throat;  that  is,  the  solution  is  singular  there  as  indi- 
cated previously  It  was  decided  to  analyze  the  flow  in  the  region 

of  the  throat  and  thus  obtain  a solution  which  is  uniformly  valid 
throughout  the  channel,  including  the  throat.  This  was  done  and  the 
results  included  in  the  solutions  used  in  making  the  movie  referred 

to  in  task  (1).  In  addition,  a paper  covering  the  derivation  of  the 

(2) 

solutions  from  which  the  movie  was  made  was  written;  it  was  pre- 
sented at  the  AIAA  16ttl  Aerospace  Sciences  Meeting,  Huntsville,  Ala- 
bama, January,  1978,  and  has  been  accepted  for  publication  in  the  AIAA 
Journal . 

As  a result  of  the  extra  work  done  in  improving  the  solutions 
and  preparing  the  publication,  it  was  not  possible  to  do  work  of  any 
significance  on  task  (3).  As  noted  in  the  description  of  task  (3), 
it  was  to  have  been  attempted  only  in  the  event  tasks  (1)  and  (2)  were 
completed  before  the  contract  ended. 


Adamson,  Jr.,  T.  C.  and  Liou,  M.  S.,  "Unsteady  Motion  of  Shock  Waves 
in  Two  Dimensional  Transonic  Channel  Flows,"  Final  Technical  Report 
prepared  for  Naval  Air  Systems  Coonand,  Report  UM014534-F,  June,  1977. 

9 

Adamson,  Jr.,  T.  C.,  Messlter,  A.  F.  and  Liou,  M.  S. , "Large  Amplitude 
Shock  Wave  Motion  in  Two  Dimensional  Transonic  Channel  Flow,"  AIAA  pa- 
per no.  78-247,  A. I. A. A.  16  Aerospace  Sciences  Meeting,  Huntsville, 
Alabama,  January,  1978. 
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(II)  RESUME  OF  CASES  SHOWN  IN  COMPUTER  MOVIE 

The  analysis  of  large  amplitude  shock  wave  motion  In  unsteady 
transonic  flow,  contained  In  reference  1,  showed  that  when  conditions 
are  such  that  oscillations  In  the  back  pressure  cause  the  shock  to 
move  up  to  the  throat  and  then  disappear  upstream,  there  are  two  pos- 
sible subsequent  motions.  In  the  first,  decreasing  back  pressure 
causes  a shock  wave  to  form  at  the  throat  and  then  move  downstream  un- 
til the  back  pressure  rises  again  forcing  the  shock  wave  to  move  back 
upstream  and  pass  through  the  throat  again;  thus,  as  this  process  re- 
peats Itself,  the  shock  wave  oscillates  about  the  throat.  In  the  sec- 
ond, decreasing  back  pressure  again  causes  the  shock  wave  to  form  at 
the  throat  and  move  downstream,  but  It  never  again  reaches  the  throat, 
oscillating  Instead  about  Its  steady  state  position.  Each  of  these 
cases  Is  Illustrated  in  the  computer  movie. 

The  movie  was  made  by  showing  the  desired  solution  at  a given 
instant  of  time  on  a screen  and  taking  a picture  of  this  solution  as 
one  frame  of  the  movie;  another  solution  at  a given  At  later  is  then 
displayed  and  another  frame  taken,  etc.  In  this  case,  the  channel 
wall,  shock  wave  position  in  the  channel,  and  the  instantaneous  pres- 
sure distribution  are  all  displayed.  The  exhaust  plenum  pressure  os- 
cillates sinusoidally  at  20  Hertz  with  an  amplitude  of  2Z  of  the  pres- 
sure at  the  throat  in  each  case.  In  order  that  the  motion  of  the  shock 
wave  and  the  corresponding  changes  in  pressure  distribution  can  be 
followed  properly,  the  solutions  are  displayed  at  1/100  real  time.  In 
each  case,  the  wall  shape  of  the  channel  is  as  follows,  where  x meas- 
ures the  distance  from  the  throat  and  is  dimensionless  with  respect  to 

i 
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the  channel  half  width. 


f(x) 


y - + (1  - 0.01  f(x)) 
w 

18x2/13  x < 1 

27(x-2)4/13  + 48(x-2)3/13  + 3 11*12 

3 x > 2 


(la) 


(lb) 


This  shape  was  the  result  of  requiring  a parabolic  form  at  the  throat 
and  continuous  first  and  second  derivatives  throughout  so  no  extrane- 
ous shock  waves  would  be  formed.  The  only  difference  between  the  two 
cases  is  the  initial,  steady  state,  location  of  the  shock  wave.  In 
case  I (oscillation  about  the  throat)  the  steady  flow  shock  wave  is 
located  such  that  the  Mach  number  entering  it  is  M » 1.08.  In  case  II, 
the  equivalent  Mach  number  is  1.10.  A brief  explanation  for  each  case 
is  contained  in  the  film. 


(Ill)  DERIVATION  OF  IMPROVED  EQUATION  FOR  CALCULATING  THE  VELOCITY 
OF  A SHOCK  WAVE  AS  IT  TRAVELS  UPSTREAM  THROUGH  A NOZZLE  THROAT 
The  solutions  for  unsteady  channel  flow  with  a shock  wave  are 
given  in  references  1 and  2.  The  source  of  the  unsteadiness  in  the 
flow  is  an  impressed  pressure  oscillation  in  the  plenum  chamber  into 
which  the  flow  exhausts;  these  oscillations  perturb  a steady  flow  in 
which  a shock  wave  occurs.  The  flow  upstream  of  the  shock  wave  is 
therefore  steady.  The  emphasis  in  references  1 and  2 is  on  large  scale 
movement  of  the  shock  wave  and  in  particular  on  those  cases  where  the 
shock  wave  actually  moves  through  the  throat,  disappears  upstream,  and 
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forms  again  at  the  throat  and  moves  downstream  when  the  back  pressure 
becomes  low  enough.  In  the  solutions  given  in  these  references,  there 
is  a singularity  at  the  throat;  for  example,  the  solutions  for  the 
shock  wave  velocity  give  infinite  shock  velocity  at  the  throat.  This 
is  clearly  physically  impossible.  As  it  turns  out  it  is  possible  to 
find  corrected  solutions  valid  in  an  inner  region  enclosing  the  throat; 
these  solutions  may  then  be  used  to  construct  uniformly  valid  composite 
solutions.  The  inner  region  and  composite  solutions  are  derived  in 
this  section. 

The  following  sketch  shows  the  channel  and  notation  used  here. 


[ 


All  lengths  are  made  dimensionless  with  respect  to  the  throat  half 

width,  £,  and  all  velocity  components  with  respect  to  the  critical 
-* 

sonic  velocity,  a . (Overbars  denote  dimensional  quantities.)  The 
slowly  varying  time  regime  is  considered  so  that  if  T is  the  time  made 
dimensionless  with  respect  to  I /a*,  then  t is  defined  such  that 

T 

T - T t T - — » 1 (la,b) 

OJZ*) 

where  T ^ is  the  characteristic  (e.g.,  the  period)  associated  with 
the  Impressed  oscillations  and 


t 


(2) 


The  wall  shape  is  arbitrary  and  given  by  (where  E « 1 here  replaces 

e in  references  1 and  2) 


y - + (1  - E2f  (x) ) (3) 

w — 

and  the  case  considered  is  that  where 

T - (k  E2)"1  (4) 

Physically,  this  case  is  that  for  which  the  time  taken  by  a signal  to 
travel  from  the  plenum  to  the  shock  wave  is  small  compared  to  the  peri- 
od of  the  impressed  oscillations,  T 

cn 

The  solutions  for  u and  u , where  u is  the  dimensionless  shock 

s s 

, (1.2) 

wave  velocity,  are 


6 


u “ 

1 + Eu^  + ■*■••• 

(5a) 

"l  ■ 

± 1 <y2i>  f<*>  + c//2 

(5b) 

U2  " 

2 * 

f"  + X,  * + h 

L K X 

(5c) 

X 

h - 

X 

1 2 C,  * G(t) 

- X [f"  + (2Y-3)u/l  + 

o 1 u. 

(5d) 

4f" 


o ,1/2  r (-13  / - — * / r \ r a/2 


c*  - -f  [(rw)c]‘'‘  z 

n«l  n 


exp.  {-mrx*/[ (y+l)Cu]  } 


cos  (nffy) 
(«  )E-1/2 

SO 


(5e) 

(5f) 


,2  d 


u ■ k E -jr  (x  + E x . + . 
s dt  so  si 


■ -E2-£gli[f[Cu  + C2d  + C<t)1 

u 


+ • • • 


(5g) 


where  C w and  are  constants  of  Integration,  with  different  values 
upstream  and  downstream  of  the  shock  wave,  and  where  the  plus  and  mi- 
nus signs  in  equation  (5b)  hold  upstream  and  downstream  of  the  shock 
wave  respectively.  The  subscripts  u and  d refer  to  conditions  imme- 
idately  upstream  and  downstream  of  the  shock  wave  respectively,  and 
is  the  value  of  u^  immediately  upstream  of  the  shock  wave  (i.e.,  eqn. 

(5b)  with  the  + sign  and  C *0);  the  term  f"  is  used  to  denote  f" (x  ). 

wo  so 

The  expansion  for  the  shock  wave  position  isx  - x + Ex  . + . . . } 

S SO  S X 

only  the  equation  for  x is  shown  and  the  value  for  x is  found  by 
integrating  equation  (5g).  Finally,  G(t)  is  proportions',  to  the  oscil- 
lation in  pressure  imposed  downstream  of  the  shock  wave  in  the  plenum; 
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it  is  therefore  an  arbitrary  function  of  time. 

The  derivation  of  the  Improved  solution  is  best  illustrated  by 

choosing  a given  wall  shape  in  the  neighborhood  of  the  throat;  here, 

the  often  used  parabolic  form  is  considered: 

f(x)  - ax2  (6) 

With  this  form,  it  is  easily  seen  that  C a X so  that  as  x ■+■  0 

J u so  so 

3 

u -*■  °°,  from  equation  (5g) , unless  the  bracket  [2yC  /3  + C0,  + G(t)] 

S U AQ  X 

(l)  30 

= C2d  + also  is  zero  at  the  given  instant.  As  it  turns  out  , 

+ G(t)  is  indeed  zero  when  the  shock  wave  is  forming  at  the  throat 
preparatory  to  moving  downstream;  analysis  shows  that  the  shock  wave 
velocity  there  is  not  singular.  However,  C^d  + G(t)  is  not  zero  when 
the  shock  wave  is  moving  upstream  through  the  throat  and  it  is  for  this 
condition  that  the  following  analysis  holds. 

When  the  shock  wave  is  very  near  the  throat,  the  fluid  velocity 
downstream  of  the  shock  wave  may  be  written  as  follows,  by  using  equa- 
tion (6)  in  equation  (5a)  to  (5f)  and  expanding  for  x <<  1: 


1 * E^~  * + E2(a(y2  - -|) 


<fr-3)  a x2  + C* 

3(y+l)  5 * 

X 


. yS 

t 2a 


(C2d  + G(t)) 


) + • • 


(7) 


Thus,  it  is  seen  that  as  the  shock  wave  moves  to  the  throat  and  so 

x 0 must  be  considered,  the  velocity  downstream  of  the  shock  also  is 

1/2 

singular  for  + G(t)  j4  0.  Moreover,  it  is  seen  when  x = 0(E  ), 

the  first  and  second  order  terms  become  of  the  same  order,  so  an  inner 

1/2 

region  must  be  considered,  of  order  E in  thickness.  In  this  region. 


— 


- ■- . 


8 


then,  a new  coordinate  x is  defined  such  that 


1/2 


x “ x/E 


and  so,  equation  (7)  may  be  written  as  follows 


(8) 


1 - e3/2  { ^ * + iW  <c2d + «*»«> 


+ E2a(,2  - i)  - E3  <2^3>  a x2  + 
E a^y  3^  E 3(y+l)  a x 


(9) 


In  the  same  way  one  can  use  equation  (5g)  to  write  the  corresponding 
equations  for  the  shock  wave  velocity: 


- *3/2  iW-  <c2d  + G<t»/I  + ••• 


(10) 


Equations  (9)  and  (10)  are  then  the  outer  solutions  to  which  the  inner 

solutions  must  match  term  by  term. 

From  equation  (9) , it  is  seen  that  in  the  inner  region 
3/2 

u - 1 + 0(E  ).  Hence,  one  can  define  a velocity  potential  where 

-1/2 

u * <!>  * E <p~  and  v ■ <t>  as  follows 

x x y 

♦ - E1/2x  + E2^x (x,y,t)  + E5/2*2(x,y,t) 

+ E3^3(x,y,t)  + ...  (11) 

Because  the  flow  is  transonic  and  the  shock  waves  weak,  it  is  possible 
to  define  a velocity  potential  to  the  desired  order. 

The  wall  boundary  conditions  are  found  by  expanding  equation  (3) 
about  x ■ 0,  writing  the  resulting  equation  in  terms  of  x,  and  requir- 
ing that  the  flow  be  tangent  to  the  walls.  Then 
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»?<*• 

± 1. 

t) 

- 0 

(12a) 

2y(i> 

± 

t) 

- + 2 ax 

(12b) 

35  <5. 

± 1. 

t) 

- 0 

(12c) 

Finally,  the  jump  conditions  are  written  relative  to  the  moving  shock 

wave^*^.  Since  from  equation  (10)  it  is  seen  that  u ■ 0(E^^)  in 

s 

the  inner  region,  the  general  expansion  for  the  shock  wave  velocity 
and  the  jump  condition  are: 

O /O 

(13a) 


u ■ E^^u  (t)  + ... 
s so 


(VUs)(ud_Us)  " 1 ' 2(£i)us  + 


(13b) 


It  should  be  noted  that  the  inner  region  in  question  is  station- 
ary, enclosing  the  nozzle  throat,  and  the  shock  wave  moves  in  it. 
Hence,  unsteady  flow  equations  hold.  For  example,  then,  the  relation- 
ship derived  for  hfc  in  the  inertial  frame  in  reference  1 is,  in  the 
present  notation. 


~2  2 2 

, a . u +v  Y+l  . / \ . . 

\ ■ * 2 ■ 2^17  + ‘VVS  + *T 


(14) 


This  is,  then,  the  equation  for  a,  the  dimensionless  speed  of  sound. 

If  equations  (8) , (11)  and  (14)  are  substituted  into  the  unsteady  flow 


JRichey,  G.  K. , and  Adamson,  Jr.,  T.  C. , "Analysis  of  Unsteady  Tran- 
sonic Channel  Flow  with  Shock  Waves,"  AIAA  Journal,  14,  1976,  pp.  1054- 
1061. 

^Chan,  J.  S.-K. , and  Adamson,  Jr.,  T.  C.  "Unsteady  Transonic  Flows  with 
Shock  Waves  in  an  Asymmetric  Channel,"  AIAA  Journal,  16,  1978,  pp.  377- 
384. 
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potential  equation linear  governing  equations  for  the  are  ob- 
tained. These  equations  are  easily  integrated  and  after  applying  the 
boundary  conditions,  equations  (12a)  - (12c),  one  finds  that  down- 
stream of  the  shock  wave,  x > x , or  x > x , 

s s 


s2  + 


i(y2 

♦l5 


(15a) 

(15b) 


where  C^(t)  and  M(t)  are  functions  of  integration.  These  solutions 
are  used  to  write  the  inner  solutions  for  u,  and  the  resulting  equa- 
tion is  written  for  x >>  1 and  compared  with  equation  (9).  For  the 
solutions  to  match  it  is  found  that  for  x > x 


C1(t)  - 2[C2d+G(t)]  (16a) 

M(t)  - 0 (16b) 


Upstream  of  the  shock  wave  the  solutions  are  similar  to  equations  (15); 
that  is,  > 0,  and  C^(t)  ■ M(t)  ■ 0 since  the  flow  is  sonic  at  the 
throat.  Hence  the  solutions  upstream  of  the  shock  wave  are  simply  the 
expansions  of  the  outer  solutions  written  in  inner  variables. 

With  known  upstream  and  downstream  of  the  shock  wave  one  can 
use  equations  (11)  and  (13)  to  show  that  the  shock  wave  velocity  in 
the  inner  region  is 


^Guderley,  K.  G. , The  Theory  of  Transonic  Flow.  Pergamon  Press, 
Addison-Wesley,  1962,  p.  7. 


(17) 


- - c&b 


u - - (■“) 

80  4 


_/wit^  + s l(t)-  i0k\ 


It  Is  seen  then  that  for  x » 1 and  as  given  in  equation  (16a),  the 
inner  shock  wave  velocity,  equations  (13a)  and  (17),  matches  with  the 
outer  shock  wave  velocity,  equation  (10),  as  it  should. 

Finally,  by  adding  the  inner  and  outer  solutions  and  subtracting 
the  common  terms  (those  used  in  matching) , one  can  write  composite 


solutions  for  u and  u for  x > x , uniformly  valid  to  order  E 
s s 


3/2 


throughout  the  channel.  Thus, 


u - 1 + E 


+ E 


{"i-  / 

> fr 

V2 


*2  + 2E<C2d  + G(t»  + 


«■} 


(y2  - b - 7(2Y-3)u?  + S\ 


♦ tcM  + c<t)l<i-+  /I I i))  + . 


(18a) 


u - - E (^) 

s 4 


{/ 


WIT  *1  * 2E(c2d  + 0(t» 


- /y+1  Xs} 


+ E 


O r>  (c2d + 5 "»  + f1 


(18b) 


It  is  seen  that  as  x -*■  0,  ug  and  u remain  finite,  as  they  should. 

It  should  be  noted  that  the  solutions  represented  by  equation 

(18)  are  necessary  only  for  C2J  + G(t)  0,  i.e.,  for  the  case  where 

the  shock  wave  is  approaching  the  throat  from  downstream,  or  when  sub- 

(1) 


As  soon  as  + G(t) 


sonic  flow  exists  throughout  the  channel, 
becomes  negative,  corresponding  to  the  shock  wave  forming  at  the  throat 
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and  moving  downstream, ^ the  shock  velocity  and  flow  downstream  of 
the  shock  are  finite  and  no  inner  solutions  are  necessary;  the  outer 
solutions  as  presented  in  reference  1 are  valid  to  the  order  Indicated. 

(IV)  SHOCK  INDUCED  SEPARATION  IN  UNSTEADY  FLOW  FIELDS 

The  fundamental  motivation  for  this  work  is  the  possibility  of 
adapting  for  use  in  unsteady  channel  or  nozzle  flows  the  known  steady 
flow  solutions  for  the  interaction  between  a shock  wave  and  a boundary 
layer. ^ In  a previous  study, ^ it  was  shown  that  the  partial  time 
derivatives  could  be  neglected  compared  to  spatial  derivatives  in  the 
conservation  equations  (e.g.,  3u/3t  <<  u9u/3x)  in  many  cases  of  tech- 

. nical  interest.  This  indicates  that  in  the  terminology  generally  used 

for  unsteady  flows,  these  problems  of  interest  fall  in  the  so  called 
slowly  varying,  or  quasi  steady,  time  regime  where  the  characteristic 
time  associated  with  the  flow  disturbance  is  large  compared  with  the 
fluid  residence  time  in  the  channel.  Since  time  derivatives  do  not 
appear  in  the  governing  equations,  the  unsteady  flow  problem  may  be 
considered  as  a series  of  steady  state  solutions  each  with  a different 
set  of  boundary  conditions.  It  should  be  noted,  however,  that  the 
proper  boundary  conditions  in  this  case  do  not  always  coincide  with 
those  used  in  solutions  to  steady  state  flow  problems;  this  point  will 
be  dealt  with  in  the  present  context  later. 

^Liou,  M.  S.,  "Asymptotic  Analysis  of  Interaction  Between  a Normal 
Shock  Wave  and  a Turbulent  Boundary  Layer  in  Transonic  Flow,"  Ph.D. 
Thesis,  University  of  Michigan,  1977. 
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In  Meiers' ^ experiments  involving  shock  induced  oscillations 
in  flow  in  a nozzle,  it  was  found  that  at  a certain  point  in  the  cycle 
a shock  occurred,  strong  enough  to  cause  the  flow  to  separate  or  to 
move  the  point  of  detachment  of  already  separated  flow.  Then,  as  the 
shock  moved  upstream,  the  separation  point  moved  with  it.  Meier  im- 
plies that  the  shock  induced  separation  is  a necessary  part  of  the 
observed  oscillatory  Instability.  However,  it  can  be  seen  in  the  ex- 
perimental pictures  that  the  distance  between  the  shock  wave  and  the 
separation  point  varies  significantly  with  time  and  grows  to  be  of  the 
order  of  the  channel  width.  If  this  behavior  is  indeed  being  induced 
by  the  shock  wave  boundary  layer  interaction  it  is  due  to  unsteady 
effects  because  such  large  distances  are  not  found  in  typical  steady 

flow  interactions  between  shock  waves  and  turbulent  boundary  layers 
(6  8) 

in  transonic  flow.  * It  was  first  believed  that  the  unsteady  effects 
were  strong  enough  to  cause  local  accelerations  to  be  important,  thus 
changing  the  problem  completely  from  that  in  the  steady  flow  case. 
However,  as  mentioned  above,  this  was  shown  not  to  be  the  case  for 
problems  of  interest  (including  Meiers'  experiments)  so  this  possibili- 
ty must  be  ruled  out.  On  the  other  hand,  it  is  clear  that  even  though 
the  solution  is  quasi  steady,  the  shock  position  and  distance  from  the 

^Meier,  G.E.A. , "Shock  Induced  Flow  Oscillations,"  AGARD-CP-168, 

Flow  Separation.  1975,  pp.  29-1  to  29-9. 

a 

Kool,  J.  W. , "Data  Report  of  a Transonic  Shock-Wave  Boundary  Layer 
Interaction  Experiment  at  a Mach  Number  of  1.44  and  1.46,"  N.L.R. 
Memorandum  AC-76-018,  The  Netherlands,  November,  1977. 


shock  to  the  separation  point  are  functions  of  time  and  it  is  of  inter- 
est to  check  to  see  if  the  predicted  values  of  these  quantities  can 
grow  to  the  size  indicated  in  the  experiments.  It  is  the  purpose  of 
this  section  to  use  known  steady  flow  solutions ^ as  a basis  for  such 
calculations  and  to  provide  an  estimate  of  the  time  for  the  formation 
and  collapse  of  a shock  wave  induced  separation  bubble.  The  results 
are  used  to  draw  conclusions  concerning  the  experimental  results. 

We  consider  a shock  wave  moving  over  a turbulent  boundary  layer 

1 

on  a flat  plate  in  transonic,  two  dimensional  flow.  The  effects  of  the 
curvature  of  the  channel,  while  certainly  important  insofar  as  quantita- 

1 

tive  results  are  concerned,  are  not  important  Insofar  as  orders  of  mag- 
nitude are  concerned  and  are  therefore  neglected.  The  flow  upstream  of 
the  shock  is  steady,  just  as  in  the  channel  flow.  We  then  take  the 
viewpoint  of  an  observer  moving  with  the  shock  wave.  The  absolute  and 
relative  points  of  view  and  the  associated  notation  are  shown  in  the 

following  sketch:  X » X + ...  is  the  shock  wave  location, 
s so 


N>  <$> 


Thus,  the  caret  is  used  to  indicate  variables  in  the  relative  frame. 
Subscripts  e,  w,  and  s refer  to  external  flow,  wall  and  shock  wave 
conditions  respectively,  while  the  subscripts  u and  d refer  to  condi- 
tions Immediately  upstream  and  downstream  of  the  shock  respectively. 

Lengths  are  made  dimensionless  with  respect  to  a characteristic  length 
I,  e.g.,  the  channel  half  width  at  the  minimum  area,  and  velocity  com- 
ponents are  made  dimensionless  with  respect  to  a*,  the  critical  sonic 
velocity  in  the  flow  external  to  the  boundary  layer.  (Overbars  denote 
dimensional  quantities.)  The  enthalpy  and  stagnation  enthalpy  are  re- 

_ 2 

ferred  to  a*  , and  the  remaining  thermodynamic  quantities  are  re- 
ferred to  their  critical  values  in  the  core  channel  flow.  The  time, 

T,  is  referred  to  I/a*. 

The  relationships  between  the  absolute  and  relative  frame  quan- 
tities are  as  follows. 

t - t (19a ,b,c,d) 

V - V (19e,f,g) 

2 

ht-h+|“  (19h,i,  j) 


where 

U_  - X + ...  (20a) 

a so 

- _ l 

T - Tt  T - Tch/ey)  (20b, c) 

a 

and  where  q and  T ^ are  the  velocity  vector  and  the  characteristic  time 
associated  with  the  impressed  flow  disturbance  respectively. 

The  flow  external  to  the  boundary  layer  is  inviscid,  with  condi- 
tions upstream  of  the  shock  wave  given.  Conditions  downstream  of  the 
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wave  are  then  found  by  application  of  the  jump  conditions  across  a 
normal  shock  wave  in  the  relative  system.  The  relative  stagnation 
enthalpy,  which  is  unchanged  across  the  shock  wave,  may  be  evaluated 
Immediately  upstream  of  the  shock  using  equations  (19i)  and  (19f)  and 
the  known  velocity  upstream  of  the  shock.  Thus 


tu 


2 2 
U _ U 

h _ n u - — m ^ - U U + — 

\ euUS  2 (Y-l)  euS  + 2 

u 


(21) 


where  since  the  flow  upstream  of  the  wave  is  steady,  h ■ (y+1)/2(y-1) » 
however,  because  the  shock  position  changes  with  time,  U u ■ Ugu(T)  and 

A 

since  U,,  * U (T) , h is  a function  of  time.  Then,  using  the  equations 

w J tU 

for  the  jump  conditions  across  a normal  shock  wave,  one  can  show  that 
the  equivalent  Prandtl  relationship  is, 

,2 

(22) 


u; 


u 


ed 


r i (u  U + — ) 1 

~ 11  (Y+l)  Vs  2 n 


U 


eu 


Mow,  for  the  case  to  be  considered  here,  conditions  for  the  chan- 
nel core  flow  coincide  with  case  (2)  of  reference  1;  that  is,  the  flow 

2 

is  transonic  and  the  shock  velocity  is  0[Ugu  - 1)  ].  The  change  in 
shock  wave  position  is  of  order  unity.  Although  the  solutions  given 
in  this  reference  are  for  the  case  where  an  oscillating  back  pressure 
causes  the  shock  wave  motion,  this  is  not  a significant  point  for  the 
present  analysis.  It  is  only  the  relative  order  of  parameters  and  the 
accompanying  form  of  the  expansions  that  are  necessary  here,  and  these 
are  independent  of  the  form  of  forcing  function  for  the  oscillation. 

In  order  to  simplify  the  notation  used  in  the  present  analysis. 
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the  flow  relative  to  the  shock  wave,  external  to  the  boundary  layer  is 

A A A 

written  as  Ugu  ■ 1 + e(T),  similar  to  the  notation  used  in  the  analysis 

of  the  interaction  between  a shock  wave  and  a boundary  layer  in  steady 
(6)  A 

flow.  The  relationship  between  e and  the  solutions  for  unsteady 
channel  flow^,  where  these  solutions  are  expanded  in  powers  of  E, 
say,  (here  E replaces  £ used  in  reference  1)  is  as  follows: 


U 


eu 


1 + E W + E \2a8o'°-T)  + ••• 


(23a) 


eu 


1 + e(T)  - o - U 
eu  S 


(23b) 


That  is,  in  general,  the  inviscid  flow  solution  for  the  channel  is  such 
that 

U(X,Y,T)  - 1 + E ux(X)  + E2u2(X,Y,T)  + ... 

and  so,  for  example,  Uu  - Uu(Xgo,Y,T);  then  equations  (21)  and  (22) 

hold  in  general  across  the  wave  with  U and  U,  replacing  U and  U .. 

u d eu  ed 

Here,  since  the  velocity  immediately  external  to  the  boundary  layer  at 

the  shock  wave  is  desired,  U is  the  value  of  U evaluated  at  the  wall 

eu 

(Y  » 0 in  the  present  notation)  as  well  as  at  X * X (T).  U . is  the 

so  sn 

2 

absolute  velocity  of  the  shock  wave  and  is  of  order  E . That  is,  for 


the  case  under  consideration  where  T * (k  E ) 


2.-1 


US  ' 


dX 
6 

dT 


so 


+ . . . 


, dx  „ 

Jr  SO  + 

R E dt 


(24a) 


dX 
s 

dt 


SO 


- 0(1) 


(24b) 


Thus  £ is  a function  of  time  defined  in  terms  of  E,  U^,  and  u^,^,... 

evaluated  at  X » X and  Y ■ 0. 
so 
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If  equation  (23b)  and  (24a)  are  substituted  into  equation  (22) 
and  only  those  terms  which  give  contributions  up  to  0(e  ) are  retained, 
then  it  is  seen  that 

V - f- u - 2<£r>y  + - «5> 

eu 

Now,  in  the  steady  flow  case^6\  the  result  is  0gd  - ujj,  with 

U ■ 1 + e.  Hence,  in  this  unsteady  case,  there  is  an  additional 
eu 

second  order  term  in  the  flow  velocity  downstream  of  the  shock  wave, 
external  to  the  boundary  layer.  Similar  terms  would  be  found  in  the 
temperature,  density,  and  pressure  in  the  external  flow. 

The  interaction  region,  as  seen  in  the  relative  coordinate  sys- 
tem,  is  shown  in  the  following  sketch.  In  the  X direction 
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its  extent  is  0(A)  while  in  the  Y direction  it  is  of  the  order  of  the 
boundary  layer  thickness ,6.  Thus,  new  variables  x and  y are  defined 
as  follows:  ^ 

A A __ 

^ ■ 5y;  y ■ Ax;  6 ■ — ; A ■ CY +l)^^e^^5  (26a,b,c,d) 

L L L 

A A A /A 

With  velocity  components  U « U(x,y)  and  V * V(x,y)  in  the  interaction 
region.  In  equation  (26),  L ■ L/l  is  the  dimensionless  length  of  the 
boundary  layer  and  is  thus  equal  to  X plus  the  distance  from  the 
channel  throat  to  the  beginning  of  the  boundary  layer;  L is  a function 
of  time.  Because  A « 1 and  is  in  fact  small  compared  to  the  extent 
of  the  inner  region  enclosing  the  shock  wave  in  the  channel  flow  solu- 
tions, the  velocity  in  the  flow  external  to  the  boundary  layer  in  the 

A A 

interaction  region  is  given  to  lowest  order  by  the  values  Ugu  and  U^, 
upstream  and  downstream  of  the  shock  respectively.  That  is,  the  inter- 
action region  is  thin  enough  that  in  this  distance  the  variations  in 
the  external  channel  flow  velocities  are  negligible. 

The  analysis  of  the  interaction  region,  in  steady  flow,  is  car- 

A (6) 

ried  out  by  considering  three  different  layers  in  the  Y direction. 

In  the  outer,  velocity  defect  region,  pressure  forces  dominate  over 
Reynolds  stresses,  so  to  the  order  considered,  id  flow  equations 

hold.  In  the  layer  immediately  adjacent  f , only  viscous  and 

Reynolds  stresses  need  be  considered  in  the  equations  of  motion  in  the 

A 

X direction;  from  the  other  equation  of  motion  one  finds  that  3P/3Y 
is  negligible,  where  P is  the  pressure.  Intermediate  to  these  two  la- 
yers is  a layer  in  which  inertia,  pressure  gradient,  and  Reynolds  shear 


stress  terms  dominate  one  equation  o£  motion  while  the  other  again 
shows  3P/3Y  negligible  to  the  order  retained.  Throughout  the  interac- 
tion region,  the  stagnation  enthalpy  remains  constant  since  the  wall  is 
taken  to  be  adiabatic  and  the  turbulent  and  laminar  Prandtl  numbers 
are  assumed  to  be  unity.  The  solutions  in  the  outer,  inviscid  flow, 
layer  may  be  found  Independently  of  the  inner  two  layers,  for  un- 
separated flow,  and  allow  calculation  of  the  wall  pressure  distribu- 
tion. With  this  pressure  distribution  known,  the  solutions  in  the  two 
inner  layers  may  be  used  to  calculate  the  wall  shear  stress  distribu- 
tion, T . The  question  now  is  whether  these  solutions  or  some  modi- 
fied form  of  them  hold  in  the  present  case. 

The  forms  of  the  conservation  equations  which  hold  in  the  rela- 
tive coordinate  system  are  derived  from  the  general  unsteady  conserva- 
tion equations  by  using  the  following  transformations 


i_  , 3_  u — -5- 

3T  M ' S U SX 


(27a, b,c) 


The  equations  which  hold  in  the  turbulent  boundary  layer  in  the  rela- 
tive system  are  then. 


+ V . pq  ■ 0 

3T 


p il  +(.oq  • V)q  + T U P V P + ( ) 

3T  sh  Y 


WU  A A A « ftp 

p + (p|  • V)h  + U U .P  • ± + ( ) 

3T  t sh  y 3T 


(28a) 


(28b) 


(28c) 


where  the  empty  brackets  in  equations  (28b)  and  (28c)  represent  viscous 
and  Reynolds  stress  terms.  Only  those  terms  which  involve  differences 
between  those  important  in  the  present  problem  and  those  important  in 

A A 

the  steady  flow  case  are  shown  here;  that  is,  since  3U/3Y  - 3U/3Y, 
Reynolds  and  viscous  stress  terms  have  the  same  form  in  either  case. 
Also,  it  should  be  noted  that  in  order  to  consider  unsteady  turbulent 
flow,  it  is  possible  to  consider  only  those  time  scales  which  are  large 
compared  to  the  time  scale  necessary  to  define  a time  average  in  the 
turbulent  flow. 

It  has  been  shown  in  reference  6 that  the  major  contributions 
to  the  expression  for  the  shear  stress  at  the  wall,  r , are  found  by 
consideration  of  the  flow  downstream  of  the  shock  wave.  This  is  due 
to  the  fact  that  for  shock  waves  strong  enough  to  cause  separation, 
the  upstream  influence  of  the  interaction  is  confined  to  a region  ex- 
ponentially small,  in  asymptotic  terms,  compared  to  the  boundary  layer 
thickness.  Hence,  it  is  sufficient  to  consider  the  undisturbed  boun- 
dary layer  flow  entering  a shock  wave  which  extends  deep  into  the 
boundary  layer,  causing  a large  variation  in  wall  pressure;  as  the 
shock  wave  is  approached  from  downstream,  a singularity  is  found.  The 
net  effect  of  this  is  that  in  the  three  aforementioned  layers  in  which 
solutions  must  be  found,  one  writes  the  solutions  in  terms  of  pertur- 
bations from  the  post  shock  conditions  in  the  external  flow.  The  same 
procedure  is  used  here. 

In  both  the  outer  invlscid  flow  layer  and  the  Reynolds  stress 
sublayer  (intermediate  layer),  then,  the  velocity  components  have 


expansions  of  the  following  form  in  the  steady  flow  analysis 


U - Ugd  + lyi^x.y)  + ...  - !-£  + ...+  ly^Cx.y)  + eiy^Cx.y) , , 


U - e1^2uTv1(x,y)  + e3^2uTv2Cx,y)  + 


with  similar  expressions  for  the  pressure,  density,  and  temperature. 


Here,  Uj,  is  the  dimensionless  friction  velocity,  where  uT  ■ 0(5),  and 

-1 


uT  « e « 1.  Again,  ■ Ueu  in  steady  flow.  The  comparable 
expression  in  the  present  case  is,  then. 


A A 


U - Ugd  + uTu1(x,y,t)  + ...  - 1 - e(t)  + ...  + Uj,fl^(x,y,t)  + 

+ uTeu2(x,y,t)  + ... 


(29) 


where  equation  (25)  has  been  used  for  Ued.  A similar  expression  holds 


for  P.  The  equivalent  expansion  for  V is 

A 

V 


e1/2uTv1(x,y,t)  + e3^2uTv2(x,y,t)  + 


(30) 


Now,  if  equations  (20b),  (26),  (29),  and  (30)  are  substituted  into  the 


X component  of  equation  (28b),  for  example,  the  result  is,  for  T 


(k  E2)”1  = [k(^~  + ...)2j_1,  and  P - 1 + £(t)  + .... 
el 
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l“T  5P- 


Now 


a/2 


, since  A » 0(u^S  ) and  6 » 0(u^,),  it  is  seen  that  to 


and 
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including  second  order  terns,  the  terms  involving  tine  derivatives  are 
small  compared  to  those  Involving  space  derivatives.  Hence,  to  the 
order  considered,  this  equation  has  the  same  form  as  that  used  in  the 
steady  flow  solution.  Similar  calculations  result  in  the  same  conclu- 
sion for  the  Y momentum  and  the  continuity  equations.  In  the  layer 
adjacent  to  the  wall  the  density  expansion  is  p ■ Pw(T)  + uTP^  + ... 

A 

with  similar  expansions  for  the  temperature  and  pressure,  and  U ■ 


u^,u^  ...  - Ug^.  Since  the  only  difference  from  the  expansions  in  the 


other  layers  is  that  U is  smaller  (of  higher  order)  the  conclusion  is 
again  that  equations  similar  to  those  used  in  the  steady  flow  case 
are  found.  The  energy  equations  gives 


tu 


-JC+1  . u + 

(Y-l)  sh  + * * * 


(32a, b) 


throughout  the  boundary  layer,  where  equations  (21)  and  (29)  have  been 

,2 


used  to  write  h to  order  § in  equation  (32b).  Again,  Uc  is  a func- 
tion  of  time  alone. 

The  Interaction  problem  becomes,  then,  a problem  with  governing 
equations  which  have  essentially  the  same  form  as  those  in  the  steady 
flow  problem.  The  only  change,  resulting  from  the  fact  that  the  stag- 
nation enthalpy  is  a function  of  time  rather  than  a constant,  is  the 


appearance  in  the  expansion  for  the  temperature  of  a term  which  is  a 


function  of  time  alone  and  0(£  ).  However,  the  boundary  conditions  are 


different  in  that  U » - Ug  at  the  wall,  and  in  the  flow  external  to 


the  boundary  layer,  the  boundary  condition  is  given  by  equations  (25) 


>-1 


rather  than  the  steady  flow  condition,  U ^ ■ (U^)  . At  each  boundary. 


then  there  is  an  added  term  in  the  velocity,  of  order  £ , which  is  a 


i 1 
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function  of  time  alone.  Because  the  wall  shear  stress  term, 

T ■ constant (fr  3U/3Y)  is  found  by  matching  solutions  for  U throughout 
w w 

the  three  layers  in  the  interaction  region ^ and  since  these  solutions 

2 

differ  from  their  steady  flow  counterparts  by  terms  of  0(g  ),  it  is 
seen  that  at  most,  the  expression  for  the  present  unsteady  flow  case 
will  differ  from  that  found  for  the  steady  flow  case  by  a term  which  is 
a function  of  time  alone  and  which  is  of  order  £ . Since  this  extra 
term  is  not  critical  insofar  as  this  analysis  is  concerned,  it  has  not 


been  derived  in  detail;  instead,  it  is  written  as 


eVo. 


Finally, 


then, by  adding  this  term  to  equation  5.79  in  reference  6,  the  steady 
flow  solution,  one  finds  the  following  equation  for  Tw: 

Tw(x)  = 1 + a t + [ a(a-l)  + F(t)]£2  + ...  uT(-  ~ P^x)) 


+ £ux(ln  A )4~  (^r  + Y + m)  + £ uT  { Pj^OO 
Co  Co  ^ 


P 

+ 2 4^  (3y  + m - ^-)(-  -^(x)) 
C C ' 


& /F 


+ 4 "~eU  a ( ——  + Y + m)  (In  x - Y + In 

C C e 

o o 


+ ... 


where 
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1 1 ^1 
Pu(x)  - (2Y-l)P1(x)  + ±(y  + y)x  ^ 


(34h) 


In  equations  (34),  Tgu  is  the  dimensionless  temperature  external  to  the 
boundary  layer  immediately  upstream  of  the  shock  wave  and  is,  there- 
fore, a function  of  time;  Y is  the  ratio  of  specific  heats;  and  uQ1 
is  the  variable  part  of  the  velocity  of  the  flow  in  the  velocity  defect 
layer  in  the  undisturbed  boundary  layer.  (U  ■ + utuq^  + •••)• 

As  shown,  for  example  in  reference  6,  a corresponding  incompressible 

..(9) 


profile  such  as  Coles 


with  the  proper  coefficients  can  be  used  here. 


The  value  for  u^  is  found  from  the  following  equation: 


C /a  In  Re  + . . . 
o 


(35a) 


Q 

Coles,  D. , "The  Law  of  the  Wake  in  the  Turbulent  Boundary  Layer," 
Journal  of  Fluid  Mechanics,  vol.  1,  1956,  pp.  191-226. 
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Re  - (Mfr  L (35b) 

0 eu 

Typical  variations  of  x^  vs  x for  the  steady  case^  (i.e.,  X8Q 

- constant,  U . ■ 0)  are  shown  in  figure  1.  This  corresponds  to  a 
sn 

plot  of  r - e2F(o  vs  x in  the  present  case,  at  a given  time  and  thus 
for  a given  e.  It  is  seen  that  x^  goes  through  a minimum.  Although 
the  solution  given  in  equation  (33)  is  an  asymptotic  expansion  and  in 
mathematical  terms  is  thus  valid  only  when  each  succeeding  term  is 
small  compared  to  the  previous  term,  it  has  been  shown ^ that  the  use 
of  this  expression  to  predict  incipient  separation  (i.e.,  ■ 0 at 

3 Tw/^x  ■ 0)  gives  results  which  appear  to  agree  with  the  limited  ex- 
perimental results  available.  Here  we  assume  that  the  solutions  give 

proper  trends  also  for  separated  flows  as  long  as  the  separation  bubble 

2 

is  small  enough.  Plots  of  x - £ F(t)  vs  x at  incipient  separation 

w 

and  when  the  flow  has  separated  to  form  a small  bubble  are  shown  in 
figures  2 and  3 respectively. ^ The  point  to  note  is  that  the  dis- 
tance from  the  shock  wave  (x  ■ 0)  to  the  point  of  separation  is  less 

than  or  equal  to  the  distance  to  the  minimum  point  of  x . Finally,  the 
2 

inclusion  of  the  i F(t)  terms  will  not  change  these  graphs  significantly, 
since  t is  at  most  0.25. 

It  is  seen,  from  figures  2 and  3,  that  as  the  Mach  number  (and 
therefore  t)  of  the  flow  entering  the  shock  wave  increases  at  constant 
Reynolds  number,  the  distance  from  the  shock  wave  to  the  separation 
point  decreases;  i.e.,  the  maximum  distance  occurs  at  incipient  separa- 
tion (figure  2).  Furthermore,  it  is  easily  shown ^ and  in  fact 
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Illustrated  In  figure  1 that  as  the  Reynolds  number  Increases  (u  de- 
creases) the  value  of  (T  ) Increases.  Thus,  for  separated  flow, 

v min 

the  slse  of  the  bubble  decreases  as  Reynolds  number  increases;  how- 
ever, it  should  be  noted  that  very  large  changes  in  Reynolds  number 
cause  only  very  small  changes  in  (Tw)min  and  thus  in  bubble  size. 

These  results  allow  the  following  interpretation  of  the  flow  in  the 
present  unsteady  channel  flow  problem.  As  the  shock  wave  moves  up- 
stream or  downstream  both  the  local  boundary  layer  thickness  upstream 
of  the  shock  (or  equivalently  the  local  Reynolds  number)  and  the  local 
Mach  number  of  the  flow  into  the  shock  wave  vary  with  time.  Since 
(Tw)min  is  a function  of  Reynolds  number  and  since  the  distance  from 
the  shock  wave  to  the  separation  point  varies  whenever  (Tw)m^n  varies, 
in  separated  flow,  it  follows  that  as  the  Reynolds  number  changes  the 
distance  from  the  shock  to  the  separation  point  changes  also.  However, 

since  (t  ) . is  relatively  insensitive  to  large  changes  in  Reynolds 
w mm 

number,  these  variations  will  be  neglected  compared  to  those  caused  by 

changes  in  Mach  number.  Furthermore,  because  the  shock  velocity  is  of 
a2 

order  e and  changes  in  the  absolute  velocity  due  to  changes  in  posi- 
tion are  of  order  £,  the  shock  wave  velocity  may  be  neglected  in 
assessing  the  most  important  effects.  Finally,  then,  the  most  important 
parameter  in  ascertaining  the  Instantaneous  distance  from  the  shock 
wave  to  the  separation  point  is  G,  which  measures  the  value  of  the 
instantaneous  absolute  Mach  number  immediately  upstream  of  the  given 
shock  wave  position. 

If  one  were  to  consider  an  unsteady  nozzle  flow  with  steady  flow 
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upstream  of  the  shock  wave  and  with  no  separation  other  than  that 

caused  by  the  shock  wave,  the  following  sequence  of  events  would  be 

predicted  by  the  above  analysis.  As  the  shock  wave  moves  downstream 

(€  increases)  separation  occurs  at  some  position  and  the  bubble  then 

proceeds  to  grow.  As  soon  as  the  shock  reverses  position  and  begins 

to  move  upstream,  the  bubble  begins  to  decrease  in  size  and  finally 

disappears.  Hence,  the  time  characteristic  of  the  foundation  and 

collapse  of  the  bubble  is  the  characteristic  time  of  the  shock  wave 

motion  Itself.  This  is  a direct  result,  again,  of  the  fact  that  so 

called  quasi-steady  solutions  hold.  The  most  important  point  to  be 

noted  here  is  that  throughout  this  sequence  of  events,  the  distance 

from  the  shock  wave  to  the  separation  point,  say  xge^ , which  is  always 

less  than  or  equal  to  the  distance  to  the  position  of  (Tw)Bjn» 

such  that  x ■ 0(1).  It  does  not  become  very  large  for  Mach  numbers 
sep 

and  Reynolds  numbers  of  interest.  Hence, 


sep 


X - X 
sep  so 


X LA*  0(A) 
sep 


(36) 


for  X * 0(1)  and  L * 0(1),  the  case  considered  here.  Finally,  then, 
so 


« 1 


(37) 


That  is,  measured  in  terms  of  the  distance  from  the  Initiation  of  the 
boundary  layer  to  the  shock  wave,  the  distance  from  the  shock  wave  to 
the  separation  point  is  small.  Although  this  result  was  deduced  from 
a solution  which  holds  only  for  small  bubbles  (at  most  of  order  6 in 
extent)  there  is  no  reason  to  suppose  that  the  first  point  of  the 
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bubble  (l.e.,  the  separation  point)  will  change  position  by  orders  of 
magnitude  because  the  tail  of  the  bubble  moves  downstream.  Indeed,  it 
appears  that  as  the  shock  wave  becomes  stronger,  the  separation  point 
moves  upstream  relative  to  the  shock  wave,  and  moves  upstream  of  the 
~ normal  shock  wave  as  a lambda  shock  is  formed;  still  the  distance  be- 
tween the  normal  shock  wave  and  the  separation  point  appears  to  be  at 
most  of  the  order  of  the  boundary  layer  thickness. 

Although  the  above  analysis  has  been  carried  out  for  the  case 
-2  2 

where  T - 0 (E  ) and  U - 0 (E  ) , it  is  seen  that  the  same  general  re- 

s 

suits  would  hold  for  the  case  t ■ 0(E  and  U ■ 0(E),  again  with 

s 

variations  of  shock  positions  of  0(1).  That  is,  shock  wave  velocities 
greater  by  an  order  of  magnitude  could  be  considered  also,  correspond- 
ing to  case  (1)  of  reference  1,  but  with  large  amplitude  shock  wave  mo- 
tions. In  this  event,  it  can  still  be  shown  that  in  the  relative  frame 
of  reference,  the  time  derivatives  may  be  neglected  to  the  desired  or- 
der, and  the  stagnation  enthalpy  is  given  by  equations  (32a, b).  The 
boundary  conditions  again  differ  from  those  used  in  the  steady  flow 

case  by  terms  involving  U • Hence,  the  result  for  t again  is  that 

s w 

given  by  the  steady  state  analysis  except  for  the  inclusion  of  a term 
involving  Us  which  is  now  of  order  e rather  than  c . The  shock  wave 
velocity  can  no  longer  be  ignored  to  lowest  order  since  it  is  0(E)  * 

A 

0(e),  but  this  leads  only  to  changes  in  numerical  values,  not  in  order. 
Evidently,  then,  the  same  conclusions  are  reached  for  this  case  as  the 

10Vidal,  R.  J.,  Wittliff,  C.  E.,  Catlin,  P.  A.  and  Sheen,  B.  H. , "Rey- 
nolds Number  Effects  on  the  Shock  Wave-Turbulent  Boundary  Layer  Inter- 
action at  Transonic  Speeds,"  AIAA  paper  no.  73-661,  AIAA  6th  Fluid  and 
Plasma  Dynamics  Conference,  Palm  Springs,  California,  July,  1973. 
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one  analyzed. 

In  summary,  che  conclusions  reached  from  the  specific  analysis 
carried  out  here  are  as  follows: 

Under  the  conditions  of  slowly  varying  transonic  flow  in  the  in- 
teraction between  a turbulent  boundary  layer  and  an  oscillating  shock 
wave: 

(a)  The  dimensionless  (with  respect  to  the  length  of  the  bounda- 
ry layer)  distance  from  the  foot  of  the  shock  wave  to  the 

1/2 

point  of  separation  is  0[(Ueu~l)  6]  where  6 is  the  order 

of  the  local  boundary  layer  thickness  (dimensionless  with 
respect  to  the  length  of  the  boundary  layer)  and  Ugu  is  the 
velocity  (dimensionless  with  respect  to  the  external  flow 
critical  velocity)  of  the  flow  external  to  the  boundary  la- 
yer, both  measured  immediately  upstream  of  the  shock  wave. 

(b)  The  time  characteristic  of  the  formation  and  collapse  of  the 
separation  bubble  is  the  same  order  as  that  associated  with 
the  oscillation  of  the  shock  wave. 

The  most  significant  difference  between  the  flow  outlined  above 
and  that  studied  experimentally  by  Meier ^ is  that  in  Meier's  flow, 
separation  was  always  present  in  the  subsonic  part  of  the  nozzle  flow, 
whether  or  not  a shock  wave  was  present.  Thus,  the  separation  bubble 
was  modified  by  but  not  caused  by  the  shock  wave.  Furthermore,  the 
separated  region  was  large  enough  that  the  resulting  core  flow  formed 
an  asymmetric  rather  than  a symmetric  channel  flow.  Nevertheless,  the 
conclusion  concerning  the  order  of  the  distance  from  the  shock  wave  to 
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the  separation  point  is  valid.  It  is  independent  of  both  the  exact 
numerical  value  of  the  pressure  downstream  of  the  thin  interaction  re** 
gion  and  the  curvature  of  the  flow  external  to  the  boundary  layer  as 
long  as  the  flow  is  transonic  as  specified. 

The  above  analysis  suggests  that  in  Meier' s^  experiments,  shock 
induced  separation  may  be  playing  a significant  part  in  the  observed 
oscillatory  instability  only  when  the  separation  point  is  very  close  to 
the  shock  wave  (within  a distance  of  order  of  the  boundary  layer  thick- 
ness just  upstream  of  the  shock  wave).  Large  scale  (of  order  of  the 
nozzle  throat  width)  oscillations  of  the  distance  between  the  shock 
wave  and  separation  point  evidently  are  not  found  in  unsteady  shock 
induced  separation  under  quasi-static  flow  conditions.  Therefore, 
under  these  conditions,  it  would  certainly  appear  that  another  mechan- 
ism, is  the  source  of  the  observed  oscillatory  instability. 

Recently,  experiments  similar  to  Meier's^,  but  more  detailed, 
have  been  carried  out  by  Sajben,  Kroutil.and  Chen^^  on  a flow  exhibi- 
ting oscillatory  instabilities.  They  identify  several  flow  regimes 
including  subsonic  flow  without  separation,  subsonic  flow  with  separa- 
tion, supercritical  flows  with  shock  waves  separated  by  large  distances 
(of  order  of  the  nozzle  throat)  from  the  separation  point,  and  super- 
critical flows  with  shock  waves  immediately  adjacent  to  the  separated 


Sajben,  M. , Kroutil,  J.  C.  and  Chen,  C.  P.,  "Unsteady  Transonic  Flow 
in  a Two  Dimensional  Diffuser,"  McDonnell  Douglas  Research  Laboratories 
Report  no.  MDRL  77-12,  1977.  (Presented  at  AGARD  Fluid  Dynamics  Panel 
Symposium  on  Unsteady  Aerodynamics,  Ottawa,  Canada,  September,  1977. 
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region.  In  the  last  mentioned  case,  the  separation  point  suddenly 
' moved  upstream  to  the  shock  foot  and  the  separation  bubble  nearly 

doubled  in  length.  The  authors  concluded  that  what  had  been  pressure 
gradient-induced  separation  in  the  subsonic  flow  had  changed  to  shock- 
induced  separation.  The  present  work  simply  adds  analytical  confirma- 
tion to  that  conclusion.  Evidently,  the  main  effect  of  the  change  is 
the  change  in  bubble  size.  It  is  interesting  to  note  that  these  au- 
thors also  concluded  that  it  should  be  possible  to  use  a quasi-steady 
analysis  to  study  the  shock  wave  boundary  layer  interaction. 

Apparently,  then,  contrary  to  the  implications  in  Meier' s^ 
work,  shock- induced  separation  is  not  a necessary  mechanism  in  the 
. formation  of  oscillatory  instabilities  in  nozzle  flows,  although  it 

may  be  an  important  factor  in  the  process.  Thus,  in  the  experimental 

I 

work  cited,  it  appears  that  when  supersonic  flow  occurs  in  the  nozzle, 
the  oscillations  are  linked  to  the  existence  of  a separated  region  of 
flow  but  this  separated  flow  bubble  is  essentially  the  result  of  the 
adverse  pressure  gradient  in  the  subsonic  part  of  the  flow.  Evidently, 
the  oscillations  of  the  shock  wave  position  and  the  bubble  size  are 
connected  through  the  unsteady  pressure  field  downstream  of  the  shock 
wave  such  that  phase  lags  between  the  shock  wave  motion  and  rate  of 
change  of  size  of  the  separated  flow  region  allow  an  undamped  oscilla- 
tion to  exist,  this  interaction  being  modified  but  not  changed  in  prin- 
ciple when  shock-induced  separation  occurs. 
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distributions  at  the  wall  in  the  outer  region  for  various  M and  Re  (from  Ref.  6) 
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